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Abstract. There is good observationnal evidence that the Steep Decay Phase (SDP) that is observed in most Swift GRBs 
is the tail of the prompt emission. The most popular model to explain the SDP is Hight Latitude Emission (HLE). Many 
models for the prompt emission give rise to HLE, like the popular internal shocks (IS) model, but some models do not, such 
as sporadic magnetic reconnection events. Knowing if the SDP is consistent with HLE would thus help distinguish between 
different prompt emission models. In order to test this, we model the prompt emission (and its tail) as the sum of independent 
pulses (and their tails). A single pulse is modeled as emission arising from an ultra-relativistic thin spherical expanding shell. 
We obtain analytic expressions for the flux in the IS model with a Band function spectrum. We find that in this framework the 
observed spectrum is also a Band function, and naturally softens with time. The decay of the SDP is initially dominated by 
the tail of the last pulse, but other pulses can dominate later. Modeling several overlapping pulses as a single broader pulse 
would overestimates the SDP flux. One should thus be careful when testing the HLE. 
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INTRODUCTION 

Most gamma-ray bursts (GRBs) observed by the Swift satellite show an early steep decay phase (SDP) in their X-ray 
light curve. It is usually a smooth spectral and temporal continuation of the GRB prompt emission, strongly suggesting 
that it is the tail of the prompt emission [1]. It is generally explained by High Latitude Emission (HLE), where at late 
times the observer still receives photons from increasingly larger angles relative to the line of sight, due to the longer 
path lenght caused by the curvature of the emitting region. These late photons have a smaller Doppler factor, which 
results in a steep decay of the flux and in a simple relation between the temporal and spectral indices a = 2 + j3, where 
F V (T) « T~ a v~P [2]. We test the consistency of HLE with the SDP by modeling the prompt emission as a sum of its 
individual pulses, including their tails. We calculate the flux for a single emission episode in the framework of internal 
shocks, and then combine several pulses to model the prompt emission. 

EMISSION OF A SINGLE PULSE 

We consider an ultra-relativistic (r> 1) thin (of width <C R/T 2 ) spherical expanding shell emitting over a range 
of radii Rq <R<R/=Rq +AR. The Lorentz factor of the emitting shell is assumed to scales as a power-law with 
radius, F 2 = T^R/Rq)'" 1 , where Tq = T(Ro). In order to calculate the flux received at any time T by the observer 
we intergrate over the Equal Arrival Time Surface (EATS; [3]), which is the locus of points from which photons 
that are emitted at a radius R, angle 9 from the line of sight and lab frame time t reach the observer at the same 
observed time T. For a shell ejected at an observer time T e j, the first photon reaches the observer at a time T e j + Tq 
with 7b = (1 +z)J?o/[2(m + l)crg]. We also define T f = T (Rf/Ro) m+1 = T Q (l +AR/R ) m+1 , whichis the last time at 
which photons emitted from the line of sight reach the observer. 

We choose for the emission spectrum the phenomenological Band function (Band et al., 1993) spectrum, which 
generally provides a good fit to the prompt GRB emission. The co-moving peak spectral luminosity is assumed to 
scale as a power-law with radius, L' v , °< (R/Rq) , where v'(R) is the peak frequency of the emitted vF v spectrum. 
Since Internal shocks is the most popular model for the prompt emission, we consider it for the following. In this 
framework, several simplifying assumptions can be made: the outflow is expected to be in the coasting phase (m = 0), 
and the electrons are expected to be in fast cooling regime. The emission mechanism is assumed to be synchrotron. 
This leads to v' p « R d with d = -1, and L' v , °c (R/R ) 1 , i.e a = 1. Then, 7b = (1 +z)R /(2cT 2 ) ), T f = 7b(l +AR/R ) 
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FIGURE 1. Left: Evolution of the shape of one normalized pulse with the normalized frequency v/vo- Middle: Evolution of 
the observed spectrum with time (corresponding to the values of T /ff written near each spectrum). The thin lines correspond to 
the rising part of the pulse, the thick lines to the decaying part of the pulse. AR/Rq = 1. v/Vo(7b) = 1. Right: Comparison of the 
evolution of the spectral (2 + 13; thin lines) and temporal (a; thick lines) slopes at fixed observed frequencies (for E'q = 0.5 keV 
and To = 300, so that E 0:Obs = 300 keV). 
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where S is the normalized Band function, x = v'/v' p , with v' = (1 +z)v/8 where v is the observed frequency, 
x b = {b\ — i»2)/(l +bi), and b\ and Z?2 are the high and low energy slopes of the spectrum; z is the redshift of the 
source and cIl the luminosity distance between the source and the observer. We define Vo = 2TqVq/(1 +z), where 
v o = v p(^o)- One should note that most of the results derived in the following hold only in the model of internal 
shocks, and not in a more general case. 

The observed flux is then (in the framework of internal shocks): 
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where^b = (1 + z)Lq / (\2ndj) . Figure 1 (left panel) shows the variation of a pulse shape with the normalized frequency 
v/vq. Different shapes can be obtained, form spiky to rouder. For m~0 and d = — 1 the observed spectrum is a pure 
Band function, just like the emitted spectrum (see middle panel of Fig. 1), where the observed peak frequency of 
the vF v spectrum decreases with time as v p = Vq/T, and T = (T — T c j)/To = 1+7'. This corresponds to a softening 
of the spectrum with time (Fig. 1 right panel) which agrees with observations. On the same panel we compare the 
evolution of the instantaneous spectral slope j3 = —d\ogF v /d\ogv with the temporal slope a = —d\ogF v /d\ogf, 
where T = (T — T e j )/Tq = 1 + 1 : we can see that the HLE relation a = 2 + j3 is valid as soon as T > Tf. One should be 
careful that this is true only in the framework of internal shocks model, and with this definition of the temporal slope 
(for exemple, a = —d\ogF v /d\ogT, which is another definition of the temporal slope, approaches 2 + J3 only at late 
time). 



COMBINING PULSES TO OBTAIN THE PROMPT EMISSION 

Within our model, the prompt emission is the sum of independent pulses, and the SDP is thus the sum of the tails of 
these pulses. For a prompt emission composed of several equal pulses, at late time the contribution of each pulse is 
equal, and the temporal slope just after the peak of a pulse increases with its ejection time T e j. When varying several 

parameters among the different pulses, the late time flux ratio of the pulse tails is the ratio of their F pea ]J^^ H . Just 
after the peak of the last pulse, the SDP is dominated by the last pulse. This shows that several pulses can dominate 
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FIGURE 2. Left: Exemple of a prompt emission consisting of three pulses with T e j = — 1 s, 13 s, 21 s, To = 2 sfor all three pulses, 
AR/Rq = 3, 2, 1, and F pea ]j/Fo = 0.7, 1, 0.7. Thin non-solid lines represent individual pulses, while the thick solid line shows the 
total prompt emission. Right: Comparison between a fit with several pulses (12 here) and a fit with one broad pulse. Thin non solid 
lines shows each individual pulse, the thin solid line shows the total prompt emission, and the thick solid line shows a possible fit 
with one broad pulse. The normalized frequency is v/vq = 0.1. Both panels are in logarithmic scale. 



the SDP at different times, as one can see in the left panel of Fig. 2. Therefore, one should be careful to consider this 
when studying the temporal and spectral behavior of the SDP. 

Figure 2 (right panel) shows what can happen if, because of noisy data or coarse time bins, a prompt emission (thin 
solid line) which is actually composed by several pulses is fitted by one broad pulse (thick solid line): the fit would 
give a tail with the same temporal slope at late time than the actual prompt tail, but with no higher temporal slopes just 
at the end of the prompt, the whole tail of this broad pulse being close to a power law. Moreover, this overestimates 
the flux of the SDP. It is important to keep this in mind when confronting such a model with actual data. 

CONCLUSION 

We have outlined a model for the prompt emission and its tail. This model contains a restricted number of free 
parameters, 10 per pulse: a, m, d, Fq, b\, b^, Eq(To), Tq, Tf and T e j. In the case for internal shocks, this can be 
reduced to 7: m = 0; d = — 1 and a = 1; as in this framework AR ~ Rq is expected, one can fix Ff/To = 1 + AR/Rq = 2. 
For a prompt emission with N pulses, the total number of free parameters can be further reduced to 3(N + 1), instead 
of 6N, as we expect the Band function parameters (b\, b^ and Eq(Tq)) to be similar for all pulses. 

The shape of a pulse can vary considerebly in our model, from very spiky to rounder, which qualitatively reproduces 
the observed diversity. The observed spectrum is a pure Band function as the emitted one in the case of internal shocks, 
and our model naturally produces a softening of the spectrum, as is observed. 

When combining several pulses to model the prompt emission, the Steep Decay Phase is initailly dominated by the 
last pulse, and is dominated at late times by the pulse with the largest Fp^Tf H (essentially the widest pulse, except 
if there is a large difference of flux between the pulses), but can be dominated by other pulses in between. 

When fitting data, one should be careful not to consider several overlapping pulses as a single broad pulse, which 
would lead to an overestimate of the prompt tail flux and a misinterpretation of the steep decay phase. 

REFERENCES 

1. O'Brien P.T., et al., ApJ 647, 1213 (2006) 

2. Kumar P. & Panaitescu A., ApJL 541, L5 1 (2000) 

3. Granot J., Cohen-Tanugi J. & do Couto e Silva E., ApJ 677, 92 (2008) 

4. Band D.,etal.,Ap7 413, 281 (1993) 



